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Abstract This article is the written version of a talk delivered at the Workshop 
on Nonlinear Dynamics and Fundamental Interactions in Tashkent and 
starts with an introduction into quantum chaos and its relationship to 
classical chaos. The Bohigas-Giannoni-Schmit conjecture is formulated 
and evaluated within random-matrix theory. In accordance to the title, 
the presentation is twofold and begins with research results on quantum 
chromodynamics and the quark-gluon plasma. We conclude with recent 
research work on the spectroscopy of baryons. Within the framework 
of a relativistic constituent quark model we investigate the excitation 
spectra of the nucleon and the delta with regard to a possible chaotic 
behavior for the cases when a hyperfine interaction of either Goldstone- 
boson-exchange or one-gluon-exchange type is added to the confinement 
interaction. Agreement with predictions from the experimental hadron 
spectrum is established. 

1. Classical and Quantum Chaos 

In order to understand in which manner classical chaos is reflected in 
quantum systems the question has been posed: Are there differences in 
the eigenvalue spectra of classically integrable and non-integrable sys- 
tems? Billiards became a preferred playground to study both the clas- 
sical and quantum case. With the arrival of computers with increasing 
power in the late seventies diagonalization of matrices with reasonable 
size became possible. The behavior of the distribution of the spacings 
between neighboring eigenvalues turned out to be a decisive signature. 
In 1979 McDonald and Kaufman performed a comparison between the 
spectra from a classically regular and a classically chaotic system [1]. As 
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Figure 1. Nearest-neighbor spacing distributions of eigenvalues for a circle (left) 
and the Bunimovich stadium (right). Taken from Ref. [1]. 



seen in Fig. 1 they observed a qualitatively different behavior between 
the nearest-neighbor spacing distribution of the circle and the stadium. 
In the first case the spacings are clearly concentrated around zero while 
they show repelling character in the second case. There were several 
authors contributing to this discussion and we mention the papers by 
Casati, Valz-Gris, and Guarneri [2], by Berry [3], by Robnik [4] and by 
Seligman, Verbaarschot, and Zirnbauer [5]. 

Very accurate results were obtained for the classically chaotic Sinai 
billiard by Bohigas, Giannoni, and Schmit (see Fig. 2) which led them 
to the important conclusion [6]: Spectra of time-reversal invariant sys- 
tems whose classical analogues are K systems show the same fluctuation 
properties as predicted by the Gaussian orthogonal ensemble (GOE) 
of random-matrix theory (RMT). K systems are most strongly mixing 
classical systems with a positive Kolmogorov entropy. The conjecture 
turned out valid also for less chaotic (ergodic) systems without time- 
reversal invariance leading to the Gaussian unitary ensemble (GUE). 

2. Random Matrix Theory 

In lack of analytical or numerical methods to obtain the spectra of 
complicated Hamiltonians, Wigner and Dyson analyzed ensembles of 
random matrices and were able to derive mathematical expressions. A 
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Figure 2. Nearest-neighbor spacing distributions of eigenvalues for the Sinai bil- 
liard with the Wigner surmise compared to the Poisson distribution. The histogram 
comprises about 1000 consecutive eigenvalues. Taken from Ref. [6] . 

Gaussian random matrix ensemble consists of square matrices with their 
matrix elements drawn from a Gaussian distribution 

p{x) = 

One distinguishes between three different types depending on space-time 
symmetry classified by the Dyson parameter (3d = 1, 2, 4 [7]. The Gaus- 
sian orthogonal ensemble (GOE, (3d = 1) holds for time-reversal invari- 
ance and rotational symmetry of the Hamiltonian 

H mn = H nrn = H nm . (2) 

When time-reversal invariance is violated and 

Hmn — [-^ ]mn > (3) 

one obtains the Gaussian unitary ensemble (GUE, (3d = 2). The Gaus- 
sian symplectic ensemble (GSE, (3d = 4) is in correspondence with time- 
reversal invariance but broken rotational symmetry of the Hamiltonian 

H^l 2 - i ]T H^a, , (4) 

7=1 

with real and symmetric and real and antisymmetric. 

The functional form of the distribution P(s) of the neighbor spacings 
s between consecutive eigenvalues for the Gaussian ensembles can be 
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approximated by 

P(S D (s) = ap D s fe exp (-bp D s 2 ) , (5) 

which is known as the Wigner surmise and reads for example in the case 
D = 2 (GUE) 

32 2 / 4 2 



P{s) = ^ s 2 exp [ s 2 ) . (6) 

If the eigenvalues of a system are completely uncorrelated one ends up 
with a Poisson distribution for their neighbor spacings 

P(s) = exp (s) . (7) 

An interpolating function between the Poisson and the Wigner distribu- 
tion is given by the Brody distribution [8] reading for the GOE case 

P(s,lo) = a(u; + l)^exp(-a^ +1 ) , a = (jj^f) , (8) 

with < oj < 1. 

Remarkably, the Wigner distribution could be observed in a number 
of systems by physical experiments and computer simulations evading 
the whole quantum world from atomic nuclei to the hydrogen atom in a 
magnetic field to the metal-insulator transition [7]. In this contribution 
we address the situation in QCD and in hadrons. 

3. Quantum Chromo dynamics 

The Lagrangian £Q CD of quantum chromodynamics (QCD) consists 
of a gluonic part £§ CD and a part £p CD from the quarks 

£QCD _ £QCD _|_ £QCD 

1 N f 

- -jF^(x)Fr(x) + J2 ™ f )tf(x) - (9) 

/=1 

with the Dirac spinor tpf, the quark mass m,f, the number of flavors Nf, 
and the generalized field strength tensor 

Fj?{x) = d»Al{x) - d v A»{x) - gf abc A%(x)A»(x) , (10) 

where the gauge field A% with the SU(3) indices a,b,c = 1, . . . , 8, the 
coupling constant g and the structure constants f a b c of SU(3) enter. The 
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main object of study is the eigenvalue spectrum of the Dirac operator of 
QCD in 4 dimensions 

$ = $ + ig^ = 7^ + i(/7^y , (H) 

with the A a the generators of the SU(3) color-group (Gell-Mann matri- 
ces) . Discretizing the Dirac operator on a lattice in Euclidean space-time 
and applying the Kogut-Susskind (staggered) prescription, leads to the 
matrix 

( m ks)""' = ^ XI [<W/i,x' r ^ ^ - $r,s'+/i r^/^ f/],^ a , (12) 
where 

^ = exp(^4£(x)y ) (13) 

are the gauge field variables on the lattice and F XI1 a representation of 
the 7 M -matrices. 

In random matrix theory (RMT), one has to distinguish several uni- 
versality classes which are determined by the symmetries of the system. 
For the case of the QCD Dirac operator, this classification was done in 
Ref. [9]. Depending on the number of colors and the representation of 
the quarks, the Dirac operator is described by one of the three chiral 
ensembles of RMT. As far as the fluctuation properties in the bulk of 
the spectrum are concerned, the predictions of the chiral ensembles are 
identical to those of the ordinary ensembles in Sect. 2 [10]. In Ref. [11], 
the Dirac matrix was studied for color-SU(2) using both Kogut-Susskind 
and Wilson fermions which correspond to the chiral symplectic (chSE) 
and orthogonal (chOE) ensemble, respectively. Here [12], we addition- 
ally study SU(3) with Kogut-Susskind fermions which corresponds to 
the chiral unitary ensemble (chUE). The RMT result for the nearest- 
neighbor spacing distribution can be expressed in terms of so-called pro- 
late spheroidal functions, see Ref. [13]. A very good approximation to 
P(s) is provided by the Wigner surmise for the unitary ensemble, 
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Pw(s) = ^s 2 e~ 4s I* . (14) 
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We generated gauge field configurations using the standard Wilson 
plaquette action for SU(3) with and without dynamical fermions in the 
Kogut-Susskind prescription. We have worked on a 6 3 x 4 lattice with 
various values of the inverse gauge coupling (3 = 6/g 2 both in the confine- 
ment and deconfinement phase. We typically produced 10 independent 
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equilibrium configurations for each (3. Because of the spectral ergodicity 
property of RMT one can replace ensemble averages by spectral averages 
if one is only interested in bulk properties. 

The Dirac operator,^* = + is anti-Hermitian so that the eigen- 
values A n of ilf) are real. Because of {^,75} = the non-zero \ n occur 
in pairs of opposite sign. All spectra were checked against the analytical 
sum rules = and X)a >o = 3V 7 , where V is the lattice vol- 

ume. To construct the nearest-neighbor spacing distribution from the 
eigenvalues, one first has to "unfold" the spectra [14]. 



Confinement (3 — 5.2 
ma = 0.05 
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Figure 3. Nearest-neighbor spacing distribution P(s) for the Dirac operator on a 
6 3 x 4 lattice in full QCD (histograms) compared with the random matrix result (solid 
lines). There are no changes in P(s) across the deconfinement phase transition. 




Figure 4- Nearest-neighbor spacing distribution P(s) for the free Dirac operator on 
a 53 x 47 x 43 x 41 lattice compared with a Poisson distribution, e~ s . 
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Figure 5. Nearest-neighbor spacing distribution P(s) for U(l) gauge theory on 
an 8 3 x 6 lattice in the confined phase (left) and in the Coulomb phase (right). 
The theoretical curves are the chUE result, Eq. (14), and the Poisson distribution, 
Pp(s) = exp(-s). 



Figure 3 compares P(s) of full QCD with Nf = 3 flavors and quark 
mass ma = 0.05 to the RMT result. In the confinement as well as in 
the deconfinement phase we observe agreement with RMT up to very 
high (not shown). The observation that P(s) is not influenced by the 
presence of dynamical quarks is expected from the results of Ref. [10], 
which apply to the case of massless quarks. Our results, and those of 
Ref. [11], indicate that massive dynamical quarks do not affect P(s) 
either. 

No signs for a transition to Poisson regularity are found. The decon- 
finement phase transition does not seem to coincide with a transition 
in the spacing distribution. For very large values of (5 far into the de- 
confinement region, the eigenvalues start to approach the degenerate 
eigenvalues of the free theory, given by A 2 = ^2j l= ism 2 (27rn fl /L fl )/a 2 , 
where a is the lattice constant, is the number of lattice sites in the 
^-direction, and = 0, . . . , — 1. In this case, the nearest-neighbor 
spacing distribution is neither Wigner nor Poisson. It is possible to 
lift the degeneracies of the free eigenvalues using an asymmetric lat- 
tice where L x , L y , etc. are relative primes and, for large lattices, the 
distribution is then Poisson, Pp(s) = e~ s , see Fig. 4. 

We have also investigated the staggered Dirac spectrum of 4d U(l) 
gauge theory which corresponds to the chUE of RMT but had not been 
studied before in this context. At [3 C « 1.01 U(l) gauge theory un- 
dergoes a phase transition between a confinement phase with mass gap 
and monopole excitations for (3 < (3 C and the Coulomb phase which ex- 
hibits a massless photon for > (3 C . As for SU(2) and SU(3) gauge 
groups, we expect the confined phase to be described by RMT, whereas 
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free fermions are known to yield the Poisson distribution (see Fig. 4). 
The question arose whether the Coulomb phase would be described by 
RMT or by the Poisson distribution [15]. The nearest-neighbor spacing 
distributions for an 8 3 x 6 lattice at (5 = 0.9 (confined phase) and at 
(3 = 1.1 (Coulomb phase), averaged over 20 independent configurations, 
are depicted in Fig. 5. Both are consistent with the chUE of RMT. 

4. Hadrons 

Taking the experimentally measured mass spectrum of hadrons up to 
2.5 GeV from the Particle Data Group, Pascalutsa [16] could show that 
the hadron level-spacing distribution is remarkably well described by the 
Wigner surmise for (3 = 1 (see Fig. 6). This indicates that the fluctuation 
properties of the hadron spectrum fall into the GOE universality class, 
and hence hadrons exhibit the quantum chaos phenomenon. One then 
should be able to describe the statistical properties of hadron spectra 
using RMT with random Hamiltonians from GOE that are characterized 
by good time-reversal and rotational symmetry. 

Combined Dala 
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Figure 6. Histograms of the nearest-neighbor mass spacing distribution for hadron 
states with same quantum numbers. Curves represent the Poisson (dashed) and 
Wigner (solid) distributions. Taken from Ref. [16]. 

In order to test this experimental finding we are comparing with the 
eigenvalues of a Hamiltonian for a realistic quark model, namely the 
Goldstone-boson-exchange (GBE) constituent quark model [17]. It in- 
cludes the kinetic energy in relativistic form 

Ho = Y. y/p?+ml (15) 
i=i 

with rrii the masses and pi the 3-momenta of the constituent quarks. 
The interaction between two constituent quarks i,j 



V(ij)=V coni (ij) + V x (ij) 



(16) 
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is given by a confinement potential in linear form 



Konf(u) = Vq + Crij 



(17) 



and a hyperfine interaction consisting of only the spin-spin part of the 
pseudoscalar-meson-exchange potentials 



Here is the distance between the quarks, Bi are the Pauli spin ma- 
trices and Aj the Gell-Mann flavor matrices of the individual quarks. 
This kind of interaction is motivated by the spontaneous breaking of 
chiral symmetry. As a consequence constituent quarks and Goldstone 
bosons should be the appropriate effective degrees of freedom at low en- 
ergies. Baryons are then assumed to be bound states of three confined 
constituent quarks with a hyperfine interaction relying on the exchange 
of the Goldstone bosons. Due to the specific flavor dependence in Eq. 
(18) a reasonable agreement between the spectra of the low lying light 
and strange baryon states calculated from the model and the experimen- 
tal spectra could be achieved. In particular the ordering of the excited 
states with respect to their parities comes out correctly as is demon- 
strated in Fig. 7. It is interesting to notice that both the experiment 
and the numerical treatment have their problems to resolve the higher 
excited states. 

In order to investigate the influence of the hyperfine interaction we also 
analyze the nearest-neighbor spacings obtained with the confinement po- 
tential without the hyperfine interaction V x and with a model consisting 
of a different kind of hyperfine interaction which is based on one-gluon 
exchange (OGE). This was traditionally used in constituent quark mod- 
els and has a flavor independent spin-spin potential. Therefore it has 
principal problems in reproducing the phenomenological ordering of the 
low lying excited nucleon states. Nevertheless, for comparison we con- 
sider here a simple version of such a model, i.e., a reparametrization of 
the Bhaduri, Cohler, and Nogami model consisting of a potential of the 
form 



V x (ij) = [Ef=i V n(rij)\[Xf + Ef=4 V K {rij)\f\J 



(18) 



V(ij) =V + C nj - 




2b 



+ 



9mimj 



s 



A 2 




(19) 



and also a relativistic kinetic energy term in its Hamiltonian [19]. The 
spectra of the low lying nucleon and delta states calculated with this 
model are inserted in Fig. 7. 
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Figure 7. Low lying nucleon (left plot) and delta (right plot) states with total spin 
and parity J™'. The left and right bars are the theoretical energies predicted from 
the GBE and OGE models as described in the text, respectively. The shaded boxes 
represent the experimental energies with their uncertainties [18]. 



In Fig. 8 we present our theoretical results of the nearest-neighbor 
spacing distribution for the nucleon and the delta. Both the hyperfine 
interaction of either Goldstone-boson exchange and one-gluon exchange 
type yield spacing distributions corresponding to the GOE. One observes 
a preference for the GOE from the linear rise at the origin while the 
other ensembles are quadratic or quartic (cf. Eq. 5). It turns out that 
the linear confinement potential alone without reproducing the spectra 
yields eigenvalues with reduced correlations between their neighbors and 
thus leading towards a Poisson distribution, as seen clearly from the 
nucleon in Fig. 8. 

5. Conclusion 

We have outlined the universal applicability of random-matrix the- 
ory and have presented own studies of quantum chromo dynamics and 
hadrons. Concerning QCD, we were able to demonstrate that the nearest- 
neighbor spacing distribution P(s) of the eigenvalues of the Dirac oper- 
ator agrees perfectly with the RMT prediction both in the confinement 
and quark-gluon plasma-phase. This means that QCD is governed by 
quantum chaos in both phases. We could show that the eigenvalues of the 
free Dirac theory yield a Poisson distribution related to regular behav- 
ior. Our investigations tell us that the critical point of the spontaneous 
breaking of chiral symmetry does not coincide with a chaos-to-order 
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Figure 8. Histograms of the nearest-neighbor spacing distribution for the nucleon 
(left plots) and the delta (right plots). The data is for Goldstone-boson exchange and 
for one-gluon exchange compared to a pure linear confinement potential of the same 
strength. Curves represent the Poisson and the GOE-Wigner distributions. 
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transition. Concerning quarks building hadrons, we employed a rela- 
tivistic quark potential model allowing for meson exchange or one-gluon 
exchange. Computing the spectrum of the nucleon and delta baryon 
indicates a spacing distribution P(s) favoring the GOE of RMT. A lin- 
ear confinement potential alone without reproducing the level ordering 
is not enough to obtain the correct fluctuations between the eigenval- 
ues. Our results are in agreement with an analysis of the experimen- 
tal mass spectrum of hadrons from the Particle Data Tables. Invoking 
the Bohigas-Giannoni-Schmit conjecture, we conclude that not only the 
quarks but also the hadrons show evidence of quantum chaos. 
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